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1 Introduction 



The influence of chaotic dynamics in spatially extended systems is a field of intensive 
research. The competition between local chaos and diffusive coupling seems to be at 
the heart of pattern formation out of a random state. A detailed knowledge of this 
mechanism seems to be necessary to understand the development of structures out of 
regular and irregular states (cf. [1] and references cited therein). These concepts link 
different fields like e.g. hydrodynamic, optical and magnetic instabilities, chemical 
reactions or biological systems. It is usually tremendous complicated to gain informa- 
tion from the basic equations of motion even if numerical simulations are concerned. 
For that reason on is forced to investigate simple model systems, a strategy which has 
proven to be fruitful in several fields of low dimensional chaotic systems [2, 3]. One 
class of suitable model systems are constituted by coupled maps which are known to 
incorporate many features of real time evolution from a phenomenological point of 
view (e.g. [4]). They are accessible not only by numerical methods with moderate 
effort but also allow for a partially rigorous treatment. Especially the weak cou- 
pling limit of hyperbolic coupled maps has been considered recently [5, 6, 7]. It has 
been shown that even in the limit of infinite system size the dynamics is dominated 
by exponentially decaying correlation functions, a property which is usually termed 
space-time-mixing. The motion is similar to the dynamics of the uncoupled system. 
Indeed it has been shown that the map lattice can be decoupled by a continuous 
transformation. Furthermore it has been suggested that pattern formation can be 
related to equilibrium phase transitions in higher dimensional spin systems using a 
thermodynamical formulation via symbolic dynamics. However no definite results are 
available at the moment. Furthermore (partially) solvable model systems are lacking 
which are suitable to study these problems. 

It is the aim of this paper to investigate a simple coupled map lattice beyond the 
weak coupling regime and to get as far as possible rigorous results. To be definite a 
lattice of maps on the circle S 1 = [0, 2tt] is considered. The circle and not the interval 
is chosen as phase space because hyperbolicity, the prerequisite to obtain rigorous 
results, is guaranteed for local expanding maps [8] . The equations of motion read 

yi+i = ( T (¥ n )) {U) = f{<Pn ] ) + tE^/^) ~ } )) |mod2vr . (1) 

L denotes the lattice size, e the strength of the coupling, and (p = (p^) E S L the 
phase space coordinates. The map / which determines the single site dynamics is 
supposed to be given by 

f(<f) — 2ip + asm(ip) |mod27r . (2) 

For different values of the parameter a G [0, 2] it interpolates between the simple case 
of the shift map (a = 0) and a strongly non-hyperbolic map which has a critical 
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point (a = 2). In the latter case it resembles strongly the structure of the Smale 
complete logistic equation. Continuity of the map T requires periodicity for the 
coupling function g. For explicit calculations the simplest choice 

g(x) = sin(x) (3) 

will be considered. A global coupling between the different lattice sites has been as- 
sumed in eq.(l). This type of coupling considerably simplifies a theoretical approach. 
It has also been suggested that such a mean field approach yields the correct de- 
scription of physically more reasonable models with short range coupling above some 
critical dimension in accordance to equilibrium statistical mechanics [4, 9]. Further- 
more globally coupled models are interesting in their own right. They may be con- 
sidered on one hand as limiting cases of models with long range interactions [10]. On 
the other hand they have been considered in the context of biology neural networks 
and Hamiltonian dynamics of coupled oscillators [11, 12, 13, 14]. Special emphasis 
has been laid on the mechanism of synchronization between the different elements. 
Analogous phenomena can be studied in the model (1). In contrast to the mentioned 
approaches the randomness is not modeled by a stochastic force but is an inherent 
property of the dynamical system. 

The plan of the paper reads as follows. Section 2 reviews the well known construc- 
tion of a symbolic dynamics for the coupled map lattice (1) on an elementary level. 
With its help a space-time-mixing stationary state can be established for sufficiently 
small coupling strength e < e e c . In the subsequent sections the two limiting cases of 
hyperbolic coupled maps (a = 0) and strongly non-hyperbolic coupled maps (a < 2) 
are discussed separately. In section 3 it is shown rigorously that the breakdown of the 
space-time-mixing regime is accompanied by a global synchronization among the lat- 
tice sites for moderate coupling strength e e c < e < e™. Furthermore space-time chaotic 
transients of tremendous length are observed which are attributed to a hyperbolic re- 
peller. The case of non-hyperbolic coupled maps is much more difficult to analyse. 
In section 4 analytical and numerical evidence is given that even for infinitesimal cou- 
pling strength complicated bifurcations arise. They cause an anomalous dependence 
of the fluctuations on the system size, an effect which has been commonly observed 
in globally coupled non-hyperbolic maps and has been termed "violation of the law 
of large numbers" [15, 16, 17, 18]. Finally comments on prospective work are given. 

2 Markov partitions for coupled map lattices 

It is the objection of this paragraph to review the well known results on the symbolic 
dynamics of system (1) on an elementary level [5]. Additionally I will set up some 
notation which will be useful in the subsequent discussion. 
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To construct the symbolic dynamics only some global features of the single site 
map / and the coupling g is needed but not the explicit expressions (2) and (3). 
Especially one demands that the single site map is monotonous. It ensures that / 
admits a Markov partition 1 [8]. For simplicity in notation let us assume a binary 
partition which is written as I = [0, 7r], h = [it, 2ir] without loss of generality. / 
maps these intervals to the whole phase space in a monotonous way, /(/,) = [0, 2tt]. 
With this partition a symbol sequence (cr ,cri, . . .), <7j G {0, 1} is assigned to every 
phase space point tp. It denotes the element I ak of the partition that contains the 
image point / fc (y?), k > 2 . For further simplification in the notation let us demand 
that the coupling obeys g(0) = and max{g'(tp)} = g'(0) = 1. Obviously the choice 
(2) and (3) shares all these properties. 

The Markov partition of the uncoupled map lattice, e = can be obtained trivially 
as a direct product Ug_ = I a (o) <S> • • • <8> I a (L-i), g_ = (a^°\ . . . , a^ L ~^). The boundaries 
of these sets are given by the " hyperplanes" tp^ = and <p^ = 7T respectively. The 
elements of this partition are mapped by T to the whole phase space in a bijective 
way. If the single site map is expanding then the partition is also generating. Any 
sequence g^,^, . . . specifies a phase space point via f] k>0 T~ k (U gLk )- The dynamics is 
equivalent to a shift operation on the two dimensional spin lattice Vq,^, . . .. 

The Markov partition obtained for the uncoupled case can be carried over to the 
case of finite coupling. The discussion is somewhat simplified if one considers the 
extended map T : R L — > IR L which develops from the original evolution equation (1) 
by suppressing the modulo operation 

M = ( T( t \\ {V) = f(r!y)\ 4- - n ( f(<rM\ - f(^) 



X n+1 



M = M } ) + t Y.9 (/(#) - ■ ( 4 ) 



Although this system has orbits which tend to infinity the original dynamics is re- 
covered via (pW = x^|mod27r. Furthermore any object in the original phase space 
S L (e.g. the partition U&) can be identified with the corresponding object in the 
extended space M L by 2n periodic continuation. The action of eq.(l) respectively (4) 
is locally invertible for moderate coupling. Using the linearization 



3L) ) =/V°) 



(5) 

one obtains for the determinant (cf. appendix A) 



det(DT(x)) = Y[f(x^) 

u=0 



n i-tEtJ+^e n i-7^ 



1 If additionally / is locally expanding, |/'| > 1 then the partition is generating. 
2 More precisely: tp £ f]k>o f~ k (^k) yields a homomorphism between the dynamics and the shift 
of symbol sequences. 
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(6) 

where the abbreviation 7^ = g'(f(x^) — f(x^)) has been used. As 7^ < 1 and 
/' > by definition it follows that 



det(DT(aO) > lif'ix^) [(1 - e) L + e(l - e)^ 1 + \e 2 ^(l - e) L ~ 2 



u=0 



(7) 



Hence the determinant does not vanish for < e < e™ = 1 which proves local 
invert ibility. In contrast to the map (1) the extended system (4) is in addition globally 
invertible because its phase space is simply connected. This property yields the main 
difference between both formulations. 

We recall the fact that the boundaries of the Markov partition are given by the 
hyperplanes x^ = respectively x^ = ir in the uncoupled case e = 0. They are the 
preimages of the hyperplanes x^ = and which separate the extended 

phase space according to 2n periodicity. It is this property which guarantees that 
the elements of the partition Ug_ are mapped to the whole phase space S L . This 
observation is used to construct the partition for finite coupling < e < e™. To this 
end let us consider for every fixed lattice site z/ G{0,...,L — 1} the preimages of the 
hyperplanes x^ = respectively x^ = 2tc. They are determined by the equations 

= H xM ^ Q (x^) := /(*<■*>) + j Y.9 " f(* M )) 

2vr = H xM ^ UQ (x^) . (8) 

For every choice x^\ /i 7^ they yield a unique solution x^ because H' xM ^ vo {x^) > 
holds for < e < t™. Therefore they define codimensions 1 manifolds in the phase 
space and slice the space in "hypercubes". These sets are mapped by T to a full cube 
[0,27r] L in a bijective way. Fig.l shows this construction in the simple case L = 2. < Fig.l 
The cubes can be labeled uniquely by symbol sequences a because they develop 
continuously from the partition of the uncoupled case. This explicit construction 
yields a partition U z of the phase space so that T : Ug_ — > S L is invertible. Hence 
T _1 (?7 2 - 1 ) H Ug^ is not empty and the set ^k>oT' k {Ug_ k ) contains at least one point for 
any sequence cr^,^, . . .. Explicit expressions for the boundaries of the partition are 
given by eqs.(8). 

To make the partition a generating one, that means that r\k>oT~ k (U 2 _ k ) contains 
exactly one point, one has to impose some expansion property. In such a case the 
diameters of the successive preimages T~ n (Ua n ) H T~^ ra_1 - ) (?7 2 - n _ 1 ) H . . . H tend to 
zero. This property is obviously fulfilled if the map T expands any points 3 . It is 



3 more precisely: the image of any S neighbourhood of a point x contains a 6' neighbourhood of 
f{x) with 5' > 5. 
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shown easily (cf. appendix B) that this global expansion property coincides with 
local expansiveness. It states that infinitesimal neighbouring points separate in the 
course of the dynamics 

Vy \\Df(x)y\\>c\\y\\, Ol . (9) 

Expansiveness may pose additional constraints on the coupling strength e. In the un- 
coupled case this property is obviously fulfilled for expanding maps, |/'| > 1. Hence it 
persists up to some critical value e e c . In the range < e < t e c the dynamics is equivalent 
to a shift in a two dimensional spin lattice (a^,^, . . .). This property proves that the 
coupled system can be decoupled by a continuous coordinate transformation 4 . The 
dynamics in this state is mixing with respect to time as well as space translations [5]. 
It is however difficult to compute the explicit value e e c from eq.(9) without specifying 
the single site map and the coupling. So we postpone the evaluation to the next 
sections. 

Let me close this paragraph by stressing briefly the opposite case of large coupling. 
Because of the infinite range of interaction any lattice sites which take the same phase 
space value, = (f^\ will maintain this property forever. Such a region of constant 
phase may be called a domain. The most trivial domain is the globally synchronized 
solution ip^ = 4> n . The dynamics of such a solution is governed by the single site 
map 4>n+i = f{<t>n)- Its stability is determined by the linearized dynamics which in 
view of eq.(l) reads 



The solution can be written as 

d L 

where the amplitudes obey 

Sa n+1 (k) = \ k f'((j) n )5a n (k), \ k=0 = 1, A M0 = 1 - eg'(0) . (12) 

Stability requires that the coefficient on the right hand side is smaller than unity on 
the average for k ^ 0. This condition determines a critical coupling strength e s c via 

exp((ln|/'(v,)|))(l-6M0)) = l (13) 

where (• • •) denotes the long time average. Synchronization sets in for e > e s c . Ob- 
viously e e c < e s c holds. The transition from the space-time-mixing regime to the 
synchronized state will be at the center of interest of the next section. 



(1-6^(0))^ + 7</(0)£ 8<pV> 



(10) 



5^ = exp (i-^-ku) 5a n (k), < k < L - 1 (11) 



4 In generally this transformation seems to be of no practical use because of its complicated 
structure. 
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3 Hyperbolic coupled maps 



The most simple case of coupled shift maps, a = has been discussed in the literature 
from several point of views (e.g. for mathematical considerations [19], for the case 
L = 2 [20]). It allows for a complete analytical treatment for moderate coupling < 
e < e™ which incorporates the transition to the synchronized state. The discussion of 
this transition is at the heart of this section. 

First of all the condition of expansivity eq.(9) can be evaluated easily. Using 
eqs.(2), (3) and (5) for a = one obtains for arbitrary vectors y in the tangent space 

(y\Df(x)\y) = 2Y,{y {v) ) 2 -jT,^^ -^ {v) ){y (v) Y 

If 2 

> (2-2e)(y\y) . (14) 

Owing to the fact that the matrix (5) is symmetric for a = the condition (9) follows 
from eq.(14) for < e < e e c = 1/2. It implies that the dynamics is equivalent to a 
shift operation in a two dimensional spin lattice £o>£i> • • • an d is mixing with respect 
to the time evolution as well as space translations 5 . 

On the other hand the stability condition for the synchronized state (13) immedi- 
ately yields a local stable uniform solution for e s c = 1/2 < e. Hence the breakdown of 
space-time-mixing is related to a global synchronization in the system. The natural 
question arises what causes the transition from the space-time-mixing regime to the 
synchronized state and whether there exist other stable solutions. 

Before we enter the analysis of the general case let us consider first the simple situ- 
ation of two coupled shift maps. It contains the essence of the subsequent analysis on 
an elementary level. If we use the extended description (4) and introduce symmetric 
coordinates i^' = ± a;( ) the system can be decoupled 

x {+ \ - 2x {+) 

x^ = 2^-) -e sin (2x(-)) . (15) 

As the variables x^ are of interest only modulo 2ir the system (15) can be considered 

modulo 2n respectively 4n. The first equation describes the chaotic motion parallel 

to the diagonal of the phase space. The dynamics perpendicular to this direction 

is governed by the second one. The corresponding map is depicted in Fig.2. For < Fig. 2 

< e < the map possesses a slope larger than unity and yields a chaotic attracting 

set. Approaching e s c from below intermittency near the fixed point x^ = sets in. 



3 On a rigorous mathematical level the latter statement requires additional investigations [5]. 
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Above the critical value the fixed point becomes stable via a pitchfork bifurcation. 
Furthermore there remains an invariant repelling Cantor set in the region bounded 
by the two unstable fixed points. The transition is usually termed a boundary crisis 
[21]. The fixed point is globally stable but chaotic transients occur. 

The dynamics of the general system (4) can be treated in a quite similar fashion 
for e s c < e < e™. To this end let us consider a neighbourhood of the diagonal 6 



A 5 = {xe 



\x 



M _ x (p) | < 5} 



(16) 



In a first step the diameter 5 is chosen in such a way that the system (4) contracts 
the neighbourhood towards the synchronized state. For x n G Ag, we obtain 



» 



» 



Pn+l x n+l\ 



2e 



2--£cos(2 



x 



(fi) — x ( v ) — -rip) 



sin [x. 



(m) 



X 



.1 



O) (") 
J-n Xn 



u(m) _ T W| 



sin(x^ 



2cos(xW-xM)+ £ cos ( 2 4 p) 

P&v,P) 



T (p) _ T 0) N 



< 



2 - 2e 



Xn 

2 . sin(5J L-2 . . . 

— cos(o*) — 1 — cos(2o*) 

L L 



| T M _ I 



;i7) 



Choosing 5* sufficiently small the prefactor becomes smaller than unity and the map 
contracts the neighbourhood As, for e e c < e < e r c 



"C 



2 - 2e 



2 N sin(5* 
— cos(d*) — — 



L — 2 
H — cos(25*) 



< 1 



(18) 



It is obvious that the same property holds for the corresponding neighbourhood A$, 
of the diagonal in the phase space S L . In terms of Fig. 2 this set is contained in the 
interval whose final points are given by the two unstable fixed points. 

In a second step let us choose 5* in such a way that the map T locally expands 
distances on the complement of As, . For if e A^ and arbitrary vector y in the tangent 
space we obtain 



(y\DT(tp)\y) = 2E(y W ) 2 -fEcos(^ ) -^ ) )(^ 



> 2E(^/ H ) 

V 

> [2-2ecos(2^)](y|y) . 



2— cos(5*) 

Lj 



(19) 



6 Only for technical reasons the subsequent discussion is performed using the extended formula- 
tion. Let me stress again that As denotes the corresponding set in the phase space S L . 
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Hence the map is expansive on if the condition 

2-2ecos(2<5*) > 1 (20) 

is valid. This set corresponds in Fig. 2 to the region where the map has slope larger 
than 1. 

Above the critical coupling strength 1/2 = e e c = e s c < e both relations (18) and 
(20) can be satisfied by some (non-unique) value 5*. As a consequence the map (1) 
contracts on As„ towards the diagonal and is expanding on the complementary set. 
It remains to investigate the invariant set which is contained in this part of the phase 
space. To this end we use again the Markov partition Ug_ which is available for < 
e < e™ — 1. The sets Vg_ = Ug_C\ AfJ yield a partition of the set under investigation. 
Because T : U„ -> S L is invertible and T(A Sr ) C A s , holds it follows that T(K) D 
and T is invertible on its image (more precisely: T : Vg_ — > T(V^) is invertible). 
Therefore T~ n (Va n ) fl . . . Vg_ Q is not empty for any symbol lattice CTq, . . . , a n . Owing 
to the expansivity of the map T on Af t the set r\k>oT~ k {Vo_ ) contains exactly one 
point of the invariant set. The dynamics on this set is again equivalent to a shift in the 
two dimensional spin lattice g^,^, . . .. Finally let me mention that the invariant set 
contains no open neighbourhood because of the expansivity of the map (cf. appendix 
B). Therefore the dynamics beyond the critical value e s c is determined by a globally 
stable synchronized state and a repelling space-time-mixing invariant set. 

Let us now study how the invariant sets influence the time evolution of the coupled 
map lattice. Below the critical coupling strength e s c the space-time-mixing state yields 
a quite trivial dynamics. The numerical simulations show a very noisy behaviour. 
Therefore let me focus on the regime e s c < e. To visualize the time evolution of 
numerical simulations some global quantity is useful. Here the absolute value of the 
" mean field" r = | exp (if(xM))/L\ (cf. section 4) has proven to be a suitable 
quantity. Its value fluctuates in the space time mixing regime and saturates to 1 
in the synchronized state. A few typical time series are shown in Fig.3 for coupling < Fig. 3 
strength e = 0.65,0.8 and systems size L = 15,21. They have been obtained from 
a randomly chosen initial condition. The time series split into a random transient 
whose duration seems to increase with the system size and a sharp relaxation towards 
the synchronized solution. The sharpness of the transition enables us to introduce a 
well defined relaxation time 

N A (ip) = min{n G N, | T n (^) G A} (21) 

depending on the initial condition ip and some neighbourhood A of the synchronized 
state (e.g. the set (16)). Because of the fast relaxation this quantity is practically 

7 The sets V„ are not empty because contains points of the form tp^ =0, fj, 7^ vq, (p^°' ) = tt 
by construction. 
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independent of the choice of A as far as this set is contained in the attracting region 
(16). The relaxation time fluctuates strongly. Its distribution is given by 

P n (A) = (6(n,N A (tp))) (22) 

where (...) denotes an average over the distribution of initial points. The quantity 

(22) has been simulated numerically for various system lengths and coupling strength 

by taking an average over 10000 randomly chosen initial conditions. The relaxation 

time, that means the set A, has been defined by the condition that the mean field 

deviates from the final value 1 by less than 10~ 3 . Fig. 4 shows some representative < Fig. 4 

results. The distribution is in all cases Poisson like with an exponential tail. The least 

square fit between the maximum of the distribution and the first box that contains 

only one member of the ensemble is indicated also. The slope of this fit yields the 

tail of the distribution P n {A) ~ C\ exp(— A^n) and coincides with the inverse of the 

mean relaxation time. The e and L dependence of this quantity has been investigated 

from the numerical data. Fig. 5 contains the dependence on the system size for various < Fig. 5 

coupling strengths. The numerical capabilities restrict the analysis to small systems 

or coupling strength far above the critical value. However in the whole range a 

exponential dependence on the system size is clearly indicated, \ €> l — exp(—La € ). 

From this observation it is evident that large systems yield huge transient times. 

The slopes in Fig. 5 allow for an analysis of the e dependence also. The evaluation 

suggests a logarithmic dependence (cf. Fig. 6) , a e ~ — lnC 2 — aln(e — e*), with a < Fig. 6 

slope a ~ 0.4. Summarizing the numerical findings we have observed a Poisson like 

distribution of relaxation times where the inverse of the mean relaxation time obeys 

A £; i ~ Ci[C 2 (e — e s c ) a ] L . Near the transition point e s c and for large system size a huge 

space-time-mixing transient is observed. 

The numerical observations can be explained in terms of the invariant sets dis- 
covered above. First of all the Poisson like distribution of relaxation times is closely 
related to the decay rate of the space-time-mixing repeller. In order o apply the 
conventional terminology of decay rates [22] it is however necessary that the invariant 
set has vanishing Lebesgue measure. Even in the simple case considered here this 
property seems to be difficult to proof which renders calculations from first principles 
impossible. If we however assume that the repeller has measure zero then the shape 
of the distribution is an immediate consequence of a finite decay rate which coincides 
with the rate A e> £ (cf. appendix C). The decay rate can be estimated by simple ar- 
guments without referring to tedious calculations. Owing to the hyperbolic structure 
and the Markov partition of the space-time-mixing state the motion on this set can 
be viewed as a stochastic process in phase space. The probability that a point falls on 
one step of iteration within the attracting region (16) of the synchronized solution is 
given by its volume relative to the whole phase space measure. It is simply estimated 
as p ~ 2n(5*) L ~ 1 /(2ti) l ~ (^/(27r)) L . The probability that the dynamics settles on 
the synchronized state after n iteration steps is simply given by the Binomial dis- 
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tribution P n ~ (1 — p) n ~ x p ~ exp(— np) which reduces to a Poisson distribution for 
small probability p. Its decay rate p yields the exponential dependence on the system 
size. The coupling strength enters via the diameter 5* of the attracting region. From 
eq.(18) or (20) it can be immediately estimated as 5* ~ (e — e s c ) 1 ^ 2 in reasonable 
agreement with the numerical data. 

The transition from the space-time-mixing to the synchronized state at e s c is ac- 
companied by long transients which are caused by a space-time-mixing repeller with 
a small escape rate. The nature of this repeller makes it difficult to observe the tran- 
sition in numerical simulations if the coupling is varied adiabatically. Furthermore 
for system size L > 30 the transients may become longer than any reasonable obser- 
vation time and no synchronization will be observed at all. In the thermodynamic 
limit L — > oo the space-time-mixing state becomes in a certain sense stable. The 
transition from the synchronized state to the space-time-mixing regime which is ob- 
served if one decreases the coupling strength is always a sharp transition because it is 
brought about by the local properties near the diagonal. Hence the synchronization 
is accompanied by a strong hysteresis effect. 

4 Non— hyperbolic coupled maps 

Let us now focus on the opposite case a < 2. At the beginning I should mention that 
to my best knowledge there are no general results available whether the single site 
map (2) is chaotic in this case. Nevertheless slight modifications of the approaches 
used e.g. in [23] and numerical simulations indicate that the map has no stable 
periodic orbits and possesses a smooth invariant density for a < 2 which develops 
a singularity ~ |(/?|~ 2 / 3 at a = 2 8 . In addition the results described below have also 
been found for the case of the "logistic map" /(y?) = 2lp(2tt — (p)/n, < ip < tt, 
f((p) = 2n — f(2n — if), n < ip < 2n which is obviously chaotic. 

Although the Markov partition persists for non-hyperbolic coupled maps the ex- 
pansivity is lost for very weak coupling. Hence even the weak coupling regime shows 
complex structures and it seems to be impossible to give a general description. But 
it has been stressed recently [15, 16, 17, 18] that there seems to exist some general 
feature among globally coupled non-hyperbolic maps on which I will focus in this 
section. They are related to an unusual system size dependence of the fluctuations 
of global quantities and have been termed "violation of the law of large numbers". It 
has been suggested that this effect is related to some hidden coherence in the coupled 
system. Partial explanations based on numerical evidence and stochastic arguments 
have been given. But there seems to exist no approach which explicitly refers to the 
non-hyperbolicity of the system. I will try to make some approach in this direction. 

8 The subsequent considerations demand for a generating partition. Its existence is guaranteed 
c.g in the case of expanding maps, |(/^)'| > 1 for some n e N. 
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Let me first briefly review the phenomenological results of numerical simulations 
(see also [15]). Owing to the global coupling of the maps the interaction is caused via 
a mean field 

r n exp(i<t> n ) := I £exp (if(^)) (23) 
which is a global quantity. In terms of this field the dynamics of eq.(l) is written as 

<P { £i = FnifbP)), F n (<p) = y? + er n sin($„ - ip) (24) 

where the interaction function F n depends on the other lattice sites via eq.(23). As 
long as the motion is chaotic the map (24) resembles the structure of a stochastically 
forced single map. A naive view of the law of large numbers suggests that the fluc- 
tuations of the mean field, A cc := (r 2 cos 2 ($)) — (rcos($)) 2 , A ss := (r 2 sin 2 ($)) — 
(rsin($)) 2 , will decrease with the system size as ~ L^ 1 . This behaviour is found in 
the hyperbolic case a < 2. But in non-hyperbolic situations strong deviations are 
observed. Fig.7 shows a few representative results for coupling strength e = 0.2 and < Fig. 7 
a = 1.98, a = 2. The fluctuations saturate at some finite value if the system size 
is increased beyond a critical value. The behaviour is more pronounced for the fluc- 
tuation A ss and seems to increase by approaching a = 2 or increasing the coupling 
strength. Furthermore this "violation of the law of large numbers" seems to occur for 
infinitesimal coupling strength in the non-hyperbolic case a = 2. But it is difficult to 
obtain conclusive results from the numerical simulations because of limitations in the 
system size. A much more refined view is obtained if the whole distribution function 
of the mean field is considered. Fig. 8 shows some typical results which have been ob- < Fig. 8 
tained from a time series of length N = 2.5 x 10 6 and a random initial condition. For 
small system size the distribution resembles a Gaussian shape. If the critical system 
size is approached the peak broadens and a double peak structure develops if the size 
is increased further. The spacing of the two maxima is responsible for the saturation 
of the moment A ss . It is furthermore instructive to look at the time evolution of the 
mean field (23). A finite part of the time series corresponding to parameter values 
used in Fig. 8 are shown in Fig. 9 . Whereas the evolution behaves random below the < Fig. 9 
critical system size an intermittent behaviour is clearly observed in the critical region. 
Beyond the critical length a hopping between states of positive and negative phases 
$ at very low transition rate is detected. It may exceed the available time span of 
the numerical simulation. The numerical results resemble strongly the phenomena of 
symmetry breaking chaos transitions or crisis induced intermittency [24, 25, 26] which 
are well known in the context of low dimensional chaotic dynamics. An explanation 
in this direction will be developed in the sequel. 

A partial analytical approach may be developed by studying the evolution of the 
full phase space distribution of the coordinates ipM which is governed by the Ruelle- 
Frobenius-Perron equation. This full description seems to be to complicated to handle 
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with. But one may resort to a formulation introduced in ref.[27]. The global coupling 
enables us to reduce the full dynamics to the evaluation of a single site map (24) 
which depends on the mean field (23). Hence the motion can be analysed completely 
by considering the probability distribution that a lattice site assumes a value ip at 
time n. It reads 

Pn(<P) :=7EMrf L rf (25) 

^ v 

where 62^ denotes the 2n periodic extension of the ^-distribution. An exact closed 
equation of motion for this density is easily written down. On one hand the density 
determines the mean field (23) via 

r n exp(i$ n ) = J p n (ip)exp(if(ip))dip . (26) 

On the other hand eq.(24) yields for the time evolution 

Pn+lM = 7EW°/)(V M )-V] 

= /^[(FnO/)(^)-¥»]p„(^)# • (27) 

Eq.(27) resembles the structure of a Ruelle-Frobenius-Perron equation. But the de- 
pendence of the mean field, that means the map F n , on the density itself turns the 
relation into a nonlinear evolution equation. On this level eq.(27) is an exact conse- 
quence of the full dynamical system (1) as long as the densities are of the form (25). 
For large system size these densities tend in a weak sense to smooth distributions for 
typical initial conditions. This suggestive observation has been checked numerically 
(cf. Fig. 11). But also mathematical rigorous proofs are available in the hyperbolic 
case [28]. Hence eq.(27) yields the appropriate description of the dynamics in the 
limit of large system size if we consider sufficiently smooth densities 9 . It is much 
more easier to analyse than the full Ruelle-Frobenius-Perron equation of the coupled 
map lattice. 

The long time behaviour is of central interest. Therefore let us concentrate on the 
discussion of stationary states and their stability. To this end let me briefly mention 
the symmetries of the coupled map lattice. The full system (1), (2), (3) is obviously 
invariant with respect to phase inversion, W, ip^ — > 2ir — ip( u \ This property carries 
over to the mean field equation (27) which is invariant with respect to tp — > 2tt — p. 
Any solution is therefore either symmetric, p^(<p) = p^\2n — p) or has an inversion 
symmetric counterpart p^\p) = p^\2^ — <p) 7^ p^\p). In the former case the 
mean field is real, $^ = 0, 7r, whereas in the latter case it takes complex conjugate 
values, = r^~\ $W = — A real mean field does not destroy the inversion 
symmetry of the map (24) (cf. Fig. 10) so that a symmetric solution stays symmetric <Fig.lO 

9 It should be remarked that eq.(27) does not depend on L explicitly. The system size enters only 
through the structure of the distribution. A smooth distribution corresponds to the limit L — > 00 
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in the course of the dynamics. 

The stationary densities p*((p) are determined by 



(28) 



where the map depends on the field 



r*exp(i<I>*) = / p*(ip) exp(i/(v?)) dip 



(29) 



Eq.(28) is formally a Ruelle-Frobenius-Perron equation for the one dimensional map 
F* o /. The existence of smooth solutions and their computation is in general a 
difficult problem. At least for symmetric densities the results of [23] suggest that such 
a solution exists whereas numerical simulations indicate the existence of asymmetric 
solutions also for appropriate coupling strengths (cf. Fig. 11). If one presupposes 
the existence of stationary smooth solutions their structure can be analysed in the 
non-hyperbolic case a = 2. In that case the point of slope zero, <p = ir, introduces 
singularities into the density which can be computed by a perturbation expansion 
[29]. The case of symmetric densities is easy to handle. Here the point of slope zero 

is mapped onto the unstable fixed point (cf. Fig. 10), a situation which is identical 

(s) 

to the Smale complete logistic equation. The density pi develops a singularity 
p[ S \ip) ~ M~ 2 ^ 3 at = and is otherwise smooth. The case of asymmetric solutions 
is a little more difficult to handle (cf. appendix D). Here the point of slope zero has a 
non-trivial trajectory, uj n = (F* o f) n (n). The density develops singularities at these 
image points , pi^ ~ \ip — oj n \~ 2 / 3 whose strengths decrease exponentially with n 10 . 
Hence the structure may become quite complicated. These results persist mainly for 
a values near to the non-hyperbolic case, 2 - a < 1. In that case the critical slope 
is non-zero but small and causes sharp peaks in the density in contrast to actual 
singularities. The symmetric solution is smooth and has a sharp maximum at the 
unstable fixed point whereas the asymmetric density develops peaks on the orbit of the 
critical point. For small coupling strength e this orbit yields a sequence which departs 
exponentially from the unstable fixed point (cf. Fig. 11 for comparison with numerical 
simulations). Finally it should be mentioned that the described enhancement of the 
density reflects a partial synchronization among lattice sites in the full system. This 
effect is attributed to the non-hyperbolicity of the single site map. 

Let us give some qualitative estimates on the stability of the stationary density. 
As long as linear stability analysis is valid it is determined by the eigenvalue problem 
which emerges from the linearization of the evolution equation (27). It reads 




The result is rigorous if the orbit terminates at some unstable periodic orbit [23]. 
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-elm (j 6'[<p - (-F* o f)^)]pM exp[-i/(^)] # 

• | exp[i/(^)]/i A (V) #) • (30) 

In order that this expression makes sense the density should be differentiable. For 
that reason let us restrict the subsequent discussion to the case of symmetric densities 
and a < 2. Integrating the eigenvalue equation (30) with respect to ip one obtains 
immediately that eigenvalues A ^ 1 are necessarily related to eigenfunctions with 
vanishing average value, / h\{ip) dip = 0. Furthermore there exists an eigenvalue A = 1 
which is independent of the coupling strength e. This "Goldstone mode" is related to 
a continuous symmetry of the fixed point equation (28) 11 . For small coupling eq.(30) 
yields a perturbation of the Ruelle-Frobenius-Perron equation of the symmetric map 
F* o f. The unperturbed eigenvalue problem possesses by presupposition (cf. the 
remarks at the beginning of this section) one eigenvalue A = 1. The remaining part of 
the spectrum has modulus smaller than 1 and is separated from the largest eigenvalue 
by a finite gap. This gap shrinks to zero if a approaches 2. The largest eigenvalue 
A = 1 persists if we turn on a small coupling e > 0. We need a finite value in order 
that additional eigenvalues cross the line Re(A) = 1 and induce the instability of the 
symmetric density. The strength of this critical coupling is supposed to tend to zero 
in the limit a j 2 because the gap vanishes in that limit. Hence in the non-hyperbolic 
case even an infinitesimal coupling e may cause the instability of the symmetric density 
because the leading eigenvalue A = 1 is degenerated with a continuous part of the 
spectrum [29, 30]. Summarizing this reasoning one may expect that the symmetric 
density, = 0, is stable up to a critical coupling strength. This critical value tends 
to zero in the limit a f 2. 

Further analysis is possible if we look at numerical solutions of the mean field 
equation (27). This analysis captures also dynamical properties of the system. To 
implement the numerical algorithm a sufficiently fine partition of the phase space 
[0, 27r] in intervals of equal size is considered. The densities p n are approximated by 
step functions on this partition. Eq.(27) then turns into a matrix equation which 
can be iterated easily. For our purpose a number of ~ 10 4 boxes is sufficient to give 
accurate results for the densities too. First of all it has been checked that the mean 
field equation yields the correct description in the limit of large system size. As a typ- 
ical example Fig. 11 displays the stationary densities for to parameter combinations <Fig.ll 
e = 0.05, a = 1.99, 1.996 in comparison with the direct numerical simulation of a huge 
coupled map lattice L = 5 x 10 5 . These densities have been obtained from the time 
evolution of the map lattice (1) and the mean field equation (27) which tend to a time 
independent state for the chosen parameter values after 100 time steps. The first set 
of parameters yields a symmetric density with a sharp maximum at the unstable fixed 

11 Denoting by p*(ip,e) a (non-normalized) solution of eq.(28) then otp*(ip, e) = p*(<p,e/a), a > 
holds. 
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point whereas for the second combination of parameters a non-symmetric density is 
observed. Its maxima escape from the unstable fixed point with an exponential rate 
as predicted by the analytical analysis. The mean field result and the direct simula- 
tion coincide within the statistical errors which are caused by finite size effects. Hence 
the system shows the two types of stationary solutions described above. Beyond the 
stationary properties the mean field approach reproduces the dynamical features also. 
Fig. 11c shows a typical time evolution of the phase of the mean field $ n obtained from 
a random initial condition, that means a constant initial density. Again the mean 
field description and the simulation coincide. Let us now study the stability of the 
symmetric solution in the weak coupling case. Because analytical approaches seem to 
be difficult to apply the results of a numerical integration of the mean field equation 
(27) are presented. For sufficiently small coupling e one finds a stable symmetric sta- 
tionary density for a < 2. If the coupling is increased this solution becomes unstable. 
Fig. 12 shows the time evolution of the phase of the mean field, $ n , that emerges from <Fig.l2 
a slightly asymmetric initial condition. If we disregard a short oscillatory transient 
then an exponential relaxation towards the symmetric solution is observed below a 
critical coupling strength for a < 2. An exponential relaxation towards an asymmet- 
ric stationary state is observed if we cross the critical coupling strength. The critical 
value decreases if one approaches a = 2. In the non-hyperbolic case a = 2 for every 
coupling strength only a stable asymmetric state (cf. Fig.l2d) was found. The hole 
scenario is consistent with the qualitative linear stability analysis given above. 

The numerical treatment of the mean field equation enables us to investigate 
regions of larger coupling strength also which seem to be inaccessible by means of an 
analytical approach. The time evolution of the system may become quite complicated. 
Typically I have observed a transition from stationary to oscillatory and time chaotic 
behaviour. As a typical example Fig. 13 contains the time series of the phase of the <Fig.l3 
mean field for increasing coupling strength and two parameter values a = 1.99, 1.95. 
Although this scenario which resembles on a phenomenological level the transition 
to chaos in low dimensional dynamical systems [3] has been observed for almost all 
combinations of parameter values it is not clear to which extent it depends on the 
details of the system under consideration. 

5 Conclusion 

For a certain class of globally coupled map lattices the explicit construction of a 
Markov partition has lead us to a detailed insight into the dynamics. For coupled shift 
maps, a = 0, all invariant sets have been discovered. Beyond the weak coupling regime 
which is governed by a space-time-mixing stationary state a global synchronization 
sets in. It is accompanied by transients whose length increase exponentially with 
the system size. The scenario extends to perturbations of the shift map, a > 0. 
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But then there might occur a gap between the onset of synchronization at e s c and 
the loss of space-time-mixing, that means the loss of local expansivity, at e e c . It is 
however difficult to observe such phenomena in numerical simulations because of the 
tremendous transients. If the coupling is increased to very large values, e ^> e™, then 
the synchronized state undergoes further bifurcations and may be destroyed again. 
But these effects seem to depend strongly on the nature of the coupling in contrast 
to the above mentioned properties. 

A different behaviour is observed for strongly non-hyperbolic coupled maps, a < 2. 
For small coupling a symmetry breaking is observed in the thermodynamic limit of 
infinite system size. This bifurcation is triggered by the non-hyperbolic points and 
can be understood qualitatively in terms of a mean field transfer operator. Fluctu- 
ations which are introduced by finite size effects cause jumps between the different 
asymmetric states. They are responsible for the bimodal structure of the distribution 
of the mean field and the anomalous behaviour of the mean square deviations. If the 
size of the system is increased the time scale of the intermittent hopping may ex- 
ceed the observation time so that the computed mean square deviation drops again. 
But then secondary bifurcations to time dependent stationary states lead to a lower 
bound for the mean square deviations. Whether the symmetry breaking of the in- 
version symmetry is a necessary condition for the subsequent bifurcations like in low 
dimensional systems [31] is an unsolved problem. Map lattices without inversion 
symmetry like the frequently discussed coupled logistic maps are believed to undergo 
these "secondary" bifurcations. This scenario will explain the so called "violation of 
the law of large numbers" . Results will be published elsewhere. 

Depending on the non-hyperbolicity of the single site map the breaking of the 
inversion symmetry precedes the appearance of a synchronized state. For intermediate 
values of the parameter a both bifurcations may interact. The consequences of such 
higher order codimension bifurcation is at the moment not clear. Furthermore it is an 
open question whether it influences the dynamics on a time scale that is observable 
in numerical simulations of large coupled map lattices. 

The model system proposed in this publication constitutes an example which can 
be handled to a great extent by analytical methods. It seems therefore promising 
to construct explicitly the two dimensional spin Hamiltonian and to analyse how the 
dynamical synchronization is related to phase transitions in the corresponding spin 
lattice. Such a relation might be helpful to understand pattern formation in more 
complicated and realistic systems. Work in this direction is in progress. 
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Note added in proof 

After the submission of the manuscript I got knowledge of ref. [32] which treats a 
related subject. I also thank the referee for pointing to this publication. 

Appendix A 

The determinant of the expression (5) reads 

L-l 

det(DT(x) = [] f(x^)det(A) (31) 

is=0 

where the matrix on the right hand side is given by 

Kii = &vn - y ^2 r ) vp b Vii + . (32) 

L p L 

For simplicity the abbreviation 

lu, ■= g' (/(*<">) - /(iM)) = cos ((/(*<">) - /(iM)) (33) 

has been used. We calculate the matrix of the expression (32) via its eigenvalue 
equation. It is written as 

Auv^ = j [c v £ + s u J2 V 00 ) (34) 

where use has been made of the abbreviations 

K : = A -l + y^T^ 
L » 

C + isM := exp(i/(a> } )) (35) 
and the trigonometric identity 

7vn = c„c p + s u s^ . (36) 

A denotes the eigenvalue and the components of the eigenvector. Introducing the 
two real quantities a, (5 

a + i(3:=J2(cn + i^)vM (37) 
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the eigenvalue equation (34) reduces to 

e ^ c v a + s v f3 



a 



C X ^ I 

= t2>" 



V 



A, 



P = Z X>^— • (38) 

The condition for a non-trivial solution of this homogeneous linear equation yields a 
polynomial of degree L in the eigenvalue A 

Yl^-jT, n A, + i(f) 2 E(l-7i) n A, = . (39) 

Eq.(39) is the characteristic polynomial of the matrix (32). Hence its value at A = 
yields the determinant of the matrix. Combining this expression with eq.(31) leads 
to the result (6). 



Appendix B 



Presuppose that the map T : R L — > IR L is sufficiently smooth and locally expansive, 
that means that eq. (9) holds. Consider an arbitrary phase space point x and let 



U s (x) := {y E 



\y 



x\ 



<5} 



(40) 



denote a 5 neighbourhood of this point. The image T(Us(x)) contains a 5' neighbour- 
hood of x' = T(x), U$'(x?). The largest neighbourhood is obtained if we choose the 
distance between the point x' and the boundary of the image as 5'. It will be shown 
that 5' > 5 holds. Denote by y' a common point of the boundaries of Us'(x') and 
T(Us(x)). Let s'(t) = tx' + (1 — t)y', t e [0, 1] be the straight line which connects 
these points. It possesses a unique preimage s{t) with s(0) = x and s(l) = y. By 
construction the final point y is contained in the boundary of the set (40). We now 
have 



5' = 



\y -£ 



dT(s) 



dt 



dt > 



Df(s(t)) 



ds 
~dl 



dt > c5 



(41) 



The last inequality follows from the condition of local expansivity (cf. eq(9)) and the 
fact that the path s(t) has at least the length 5. 

Denote by T the extension of the coupled map lattice and by A the 2ir periodic 
continuation of the invariant set on which the map is locally expansive. Suppose A 
contains some neighbourhood. We will show that it coincides with the whole phase 
space. To this end suppose that for some point x G A there exists a 5 > so that 
Us(x) C A holds. Because of the invariance of A every image of the neighbourhood 
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is contained in the invariant set, T n (Us(x)) C A. Hence we can apply the global 
expansion property which states that the set T n (Us(x)) contains a neighbourhood 
with a diameter c n 5. Choosing n large enough it can be managed this set contains 
a full hypercube [0,27r] L . Hence the invariant set A contains a full hypercube. But 
then this 2tt periodic set coincides with the whole phase space. 



Appendix C 

Let A denote a neighbourhood of the attractor. Then eqs.(21) and (22) immediately 
lead to 

n 

J2Pk(A)=\(T- n (A)) (42) 

fc=0 

where A denotes the normalized measure which yields the average in eq.(22). A natural 
choice for this measure is the normalized Lebesgue measure which corresponds to a 
uniform distribution of initial conditions. Let us presuppose that the quantity (42) 
tends to 1 in the limit n — > oo. It means that the system has only one attracting set 
and that all repellers have (Lebesgue) measure zero. From the normalization of the 
measure and the identity T~ n (S L / A) = S L /T~ n (A) one obtains 

oo 

]T Pk(A) = \(T~ n (S L /A)) . (43) 

fc=n+l 

The conventional definition of the escape rate from a repelling set A is given by [22] 

a = - lim lim - In \(T' n (U e ) (44) 

where U e denotes an e neighbourhood of the repelling set. It is assumed that the escape 
rate does not depend on the choice of the (sufficiently small) neighbourhood so that 
one limit can be suppressed. The set S L / A constitutes a suitable neighbourhood. 
Then eq.(44) yields the asymptotic behaviour 

\{T- n {S L /A)) ~C ie xp(-n(7) . (45) 

Eqs.(43) and (45) result in the asymptotic relation 

P n (A)~C 2 exp(-na) . (46) 



Appendix D 

Eq.(28) can be viewed as a Ruelle-Frobenius-Perron equation for the asymmetric 
map F* o /. It may be solved by iterating a smooth initial density h n= o((p) until a 
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stationary state is reached [29] 

K+M = J S 2 A(F^ o f)ty) - V }h n (ij) dij . (47) 

For later use the explicit dependence of the map on the phase of the mean field, 
$*, has been indicated. Evaluating the integral the relation can be written as 

OV/S 4 h n [2n-(F^of)-\2n-ip)] 
n+im o o /)-i(v,)] + o o /)-i(2tt - </>)] 1 J 

where the the function / is restricted to the interval [0, tt]. If the initial function 
h n= o((p), is smooth the first iterate develops singularities at the critical points of the 
denominator. They are determined by 

(F^of)-\ v )U & (flu! 
(F^of)-\27r-ip)U & xlu! . (49) 

where io\ := er* sin(<I>*) denotes the first image of the critical point. The strength of 
the singularity can be evaluated easily by considering neighbourhood of this point. 

i) <f I uj\\ Straightforward Taylor series expansion yields to the leading order 

(^o = Q + 0(2) 

(^•.o/)'[(F..o/)- 1 ( ¥ ,)] = 4^(0) + 0(1) 

(F_^o/)- 1 (27r-^ = ^-(s^-^y) +0(2/3) 

/ - \ 2/3 

(F^of)>[(F^of)-\2n- V )] = F ^M[3jr^) +0(1) (50) 

where O(n) denotes contributions of order (ip — uoi) n . Then eq.(48) reads 

>[(y-^)/«(0))] MvO fl , 

n+AV) ~ 4^,(0) + F^(0)[3( V ,-c 1 )/^(0p/3 • 

ii) ip t uji. An analogous expansion leads to 

/ - \ 1/3 
(F^of)-^) = vr- 3^_£ +0(2/3) 



n.(o). 



(^o/y^o/)- 1 ^)] = ^(o) +0(1) 

(F_^o / )- 1 (2vr-v.) = ^7^ + 0(2) 
(F_^o/)'[(F_^o/)- 1 (27r-^] = 4F' (0) + O(l) (52) 



2/3 
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and 

. r x W , ^[2^-^-^/(4^(0))] 

n n+lM - Fkm(iUi _ ip)/nM 2 / 3+ 4*%,(0) • W 

Eqs.(51) and (53) clearly show that the first iterate develops a singularity \tp — 
uji\~ 2 ^. On further iteration this singularity is carried along the trajectory of uj\. 
With 0J2 = o f)(uji) eq.(47) yields in the vicinity of the image point 

K +1 (<p) ~ |5 27r [(F,o/)(^)-^]|^-^|- 2 / 3 # 

~ |(F,o/)Vi)|- 1/3 |^-^r 2/3 • (54) 

Hence the strength of the singularity decreases by a factor |(F* o /)'(u;i)| -1 / 3 ~ 4 _1//3 
along the trajectory uj n . 
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Figure captions 



Fig.l Diagrammatic view of the Markov partition U a (o a (i) for two coupled maps (L = 
2) in the extended phase space. The dotted lined indicates the partition for the 
uncoupled case e = 0, and the broken line the partition for finite coupling e > 0. 

Fig. 2 Diagrammatic view of the map /(-)(a;(-)) = 2x(-) - esin(2a;(->) | mod27r for 
e > e s c . The box indicates the region which contains the repelling Cantor set. 

Fig. 3 Time series of the amplitude of the mean field for different coupling strengths 
and system sizes, (a) e = 0.65, L = 15, (b) e = 0.65, L = 21, (c) e = 0.8, 
L = 15, and (d) e = 0.8, L = 21. 

Fig. 4 Distribution of the relaxation times for e = 0.8 and (a) L = 12, (b) L = 15. 
The broken line indicates the least square fit. It is shifted by a factor two for 
clarity. 

Fig. 5 Dependence of the inverse of the mean relaxation time on the system length for 
several coupling strengths. 

Fig. 6 Dependence of the slopes in Fig. 5 on the coupling strength e — e s c . The full line 
indicates the least square fit. 

Fig. 7 Mean square deviations A cc ( ) and A ss (____) in dependence on 

the system size for two values of the nonlinearity parameter a = 1.98 (x) and 
a = 2 (□). 

Fig. 8 Distribution function of the mean field for a = 1.98, e = 0.2 and several values 
of the system size, (a) L = 320, (b) L = 640, and (c) L = 1280. 

Fig. 9 Time evolution of the mean field corresponding to the parameter values chosen 
in Fig.8, a = 1.98, e = 0.2 and (a) L = 320, (b) L = 640, (c) L = 1280. The 
upper chain of dots indicates the real part r n cos($ n ) and the lower chain the 
imaginary part r ra sin($ n ). 

Fig. 10 Diagrammatic view of the stationary mean field map o / for (a) = and 
(b) > 0. The broken line indicates the orbit of the critical point. 

Fig. 11 (a) Symmetric distribution of the phase space coordinates for a = 1.99 and 
e = 0.05. The full line indicates the result of a simulation of a map lattice of 
size L = 5 x 10 5 whereas the broken line shows the mean field distribution. 
The latter curve is shifted by a factor 1/3 for clarity because both curves would 
coincide otherwise. 10 3 boxes (simulation) respectively 10 4 boxes (mean field) 
are used to figure the distribution. 



25 



(b) Asymmetric distribution for a = 1.996, e = 0.05 on a double logarithmic 
scale. The lower broken line indicates the mean field result whereas the upper 
broken line is shifted by a factor 3. 10 4 boxes (simulation) respectively 4 x 10 4 
boxes (mean field) are used to figure the distribution. 

(c) Time dependence of the phase of the mean field for parameter values used 
in Fig. lib, a = 1.996, e = 0.05. The full line indicates the mean field result 
whereas the broken line the simulation of the map lattice (L = 5 x 10 5 ). 

Fig. 12 Time evolution of the phase of the mean field for several values of the coupling 
strength and nonlinearity parameter, (a) a = 1.97, (b) a = 1.98, (c) a = 1.99, 
and (d) a = 2. 

Fig. 13 Time evolution of the phase of the mean field for larger coupling strength. 

(a) a = 1.99. The coupling strength is chosen as (from bottom to top) e = 0.1, 
0.125, 0.15, 0.175, 0.2. The curves are shifted against each other by an offset of 
0.1. 

(b) a = 1.95. The coupling strength is chosen as (from bottom to top) e = 0.4, 
0.5, 0.6, 0.7, 0.8. The curves are shifted against each other by an offset of 0.4. 
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